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ABSTRACT 

We present strong evidence for dual active galactic nuclei (AGN) in the z = 0.36 galaxy COSMOS 
J100043. 15+020637.2. COSMOS Hubble Space Telescope (HST) imaging of the galaxy shows a tidal 
tail, indicating that the galaxy recently underwent a merger, as well as two bright point sources 
near the galaxy's center. Both the luminosities of these sources (derived from the HST image) and 
their emission line flux ratios (derived from Keck/DEIMOS slit spectroscopy) suggest that both are 
AGN and not star-forming regions or supernovae. Observations from zCOSMOS, Sloan Digital Sky 
Survey, XMM-Newton, Very Large Array, and Spitzer fortify the evidence for AGN activity. With 
HST imaging we measure a projected spatial offset between the two AGN of 1.75 ± 0.03 hr 1 kpc, 
and with DEIMOS we measure a 150 ± 40 km s _1 line-of-sight velocity offset between the two AGN. 
Combined, these observations provide substantial evidence that COSMOS J100043. 15+020637.2 is a 
dual AGN in a merger-remnant galaxy. 

Subject headings: galaxies: active - galaxies: nuclei 



1. INTRODUCTION 

In the standard A cold dark matter paradigm of 
structure formation, more massive galaxies are assem- 
bled from smaller ones in a series of merger events. 
Nearly every galaxy hosts a central superrn assive black 
hole (SMBH) (jKormendv fc Richstonil 19951 ). which im- 
plies that a merger between two galaxies nearly al- 
ways results in a merger-remnant galaxy containing two 
SMBHs. Drag from dynamical friction causes the two 
SMBHs to inspiral toward the center of the merger- 
remnant. The SMBHs spend ~ 100 Myr at separations 
> 1 k pc (jBegelman et al.l 11980b iMilosavljevic fc Merritt] 
|2001|) , then form a parsec-scale binary and ultimately co- 
alesce into a single central SMBH in the merger-remnant 
galaxy. This final coalescence is necessary to preserve the 
tight observational correlation between the mass of the 
black hole and the velocity d ispersion, or total mass, o f 
the host galaxy stellar bulge (jFerrarese fc Merrittll2000T ). 

Although SMBH pairs are a natural consequence of 
galaxy mergers, there have been few unambiguous de- 
tections of galaxies hosting SMBH pairs. If sufficient gas 
accretes onto both SMBHs, they may each be visible as 
an active galactic nucleus (AGN). To date, there have 
been definitive detections of only four galaxies hosting 
such AGN pairs. First, radio signatures of AGN activity 
in the z = 0.055 elliptical galaxy 0402+379 s how that 
it hos ts binary SMB Hs se parated by 5 h" 1 pc (Ku et al.1 
11994 iManess et alj|2004; iRodriguez et al.ll2006[) . In ad- 
dition, X-ray detections of a dual AGN in the z — 0.024 
ultraluminous infrared galaxy NGC 6240 indicate it hosts 
two S MBHs separated by 0.5 h~ l kpc (|Komossa et al.l 
2003). Finally, optical spectroscopic signatures of dual 
AGN in the red galaxies EGSD2 J142033.6+525917 at 



z = 0.71 and EGSD2 J141550.8+520929 at z = 0.62 
show these galaxies host dual SMBHs at se parations of 
0.84 fc" 1 kpc and 1.6 h~ x kp c, respectively (|Gerke et al.l 
I2007t IComerford et al.ll2009l). A fifth p ossible example 
has been proposed by iBoroson fc Lauei] (|2009D . though 
it is likely to b e an ob ject of a diffe rent nature (e.g., 
Chomock et~aT]|2009allbl: lGaskel]||2009l lLauer fc Borosonl 
2009HWrobel fc Laorll2009l) ~ 



Here we present evidence for a 1.75 ± 0.03 hr 1 kpc 
projected spatial separation and 150 ± 40 km s _1 line- 
of-sight velocity separation SMBH pair, visible as the z = 
0.36 dual AGN system COSMOS J100043. 1 5+020637.2 
in the COSMOS field dScoville et al.ll2007bh . The can- 
didate was found serendipitously while visually inspect- 
ing postage stamp images of COSMOS galaxies with 
high Sersic indices in the ACS-GC catalog (Griffith et 
al., in preparation), which includes morphology measure- 
ments for over half a million sources from five large Hub- 
ble Space Telescope Advanced Camera for Surveys (HST 
ACS) imaging datasets. In addition, the candidate has 
been previously cla ssified as an AGN by both COSMOS 
dGabor et all l2009f) and the Sloan Digital Sk y Survey 
CSDSS'I (|York et al.ll27)00l : iRichards et al.l[2002T ). We as- 
sume a Hubble constant H = 100 h km s" 1 Mpc" 1 , 
f2 m = 0.3, and f^A = 0.7 throughout, and all distances 
are given in physical (not comoving) units. 

2. OBSERVATIONS 

2.1. COSMOS Imaging 

We originally identified COSMOS 

J100043.15+020637.2 as a dual AGN candidate 
from its HST F814W ACS image taken for COSMOS 
dScoville et al.l[2007al) . This image, shown in Figure [TJ 
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Fig. 1— HST F814W ACS image of COSMOS 
J100043.15+020637.2, where North is up and East is to the left. 
The galaxy's tidal tail strongly suggests it has recently undergone 
a merger, and the two bright nuclei near the galaxy center appear 
to both be AGN. The nuclei are separated by 0'.'497 ± 0'.'009, or 
1.75 ± 0.03 h- 1 kpc. 

shows a disturbed galaxy with a long tidal tail that 
suggests the galaxy has recently undergone a merger. 
The nucleus of the galaxy contains two bright point 
sources, and we use 0'.'25 (0.88 h -1 kpc) radius apertures 
to measure the magnitu de and luminosity of eac h source 
with Source Extractor (|Bertin &; Arnoutslll996l ). 

We measure F814W apparent magnitudes of 20.0 for 
the northern source and 20.1 for the southern source. For 
comparison, the entire object has an ap parent magnitude 
of 18.3, as measured using GALFIT (jPeng et all 120021 : 
Griffith ct al., in preparation). The magnitudes of the 
point sources correspond to luminosities of 3.1 x 10 43 h~ 2 
erg s _1 (8.0 x 10 9 h~ 2 L ) for the northern source and 
2.7 x 10 43 hr 2 erg s" 1 (7.1 x 10 9 hr 2 L e ) for the south- 
ern source. Supern ovae this luminous are extremely rare 
(|Smith et al.ll2007t ), suggesting that the sources are most 
likely AGN and not supernovae. However, the HST im- 
age was taken on UT 2004 March 16 and a more recent 
high-resolution image would resolve this question defini- 
tively. Regardless, other observations strongly support 
the interpretation that both sources are AGN (§ 12. 3[) . 

With Source Extractor we also find that the pro- 
jected separation between the barycenters of the sources 
is 0'.'497±0"009, or 1.75±0.03 hr 1 kpc, and the barycen- 
ters of the sources are aligned along a position angle 
6 = -40°. 4 East of North. 

2.2. DEIMOS Slit Spectroscopy 

We used the DEIMOS spectrograph on the Keck II 
telescope to obtain a 200 s spectrum of the object with a 
600 lines mm -1 grating at twilight on UT 2009 April 23. 
The spectrum spans the wavelength range 4730 - 9840 A, 
and the position angle of the slit was 9 = —36°. 7 East of 
North. The main purposes of the slit spectroscopy were 
to verify that both central point sources were AGN and 
to measure the spatial and velocity separations between 
the two AGN emission components. The AGN emission 
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Fig. 2. — Two-dimensional DEIMOS spectrum of 

J100043. 15+020637. 2 with night-sky emission features subtracted. 
The spectrum has been smoothed by a smoothing length of 2 pixels, 
and AGN line emission at [O III] A5007 (top), Ha (middle), and 
[N II] A6584 (bottom) is shown. In each panel, the vertical axis 
spans 3"41 (12.0 h~ 1 kpc at the z = 0.36 redshift of the galaxy) in 
spatial position along the slit and the horizontal axis spans 30 A in 
rest-frame wavelength centered on the emission feature. The color 
bars provide scales for the flux in counts hour -1 pixel -1 . Each 
emission feature has two components, likely corresponding to two 
distinct AGN, and the velocity separation between the two is 150 
±40kms -1 . 

is clearly visible in [O III] A5007, Ha, and [N II] A6584, 
as shown in Figure [5] We determine the projected spa- 
tial separation for each of these three emission features 
by measuring the spatial centroid of each emission com- 
ponent individually. 

To measure the spatial centroid of an emission com- 
ponent, first we center a 2 A (rest- frame) wide window 
on the emission. At each spatial position we then sum 
the flux, weighted by the inverse variance, over all wave- 
lengths within the window. We use the spatial position 
where the summed flux is maximum as the center of a 
10 pixel window (where 1 DEIMOS pixel spans 0'/ 11) 
around the emission. We then fit a quadratic to the 
summed flux to locate a peak, define a narrow window 
centered on the peak flux, and finally compute the line 
centroid within this window. We derive the error on this 
spatial centroid by repeatedly adding noise to the spec- 
trum drawn from a Gaussian with variance matching the 
DEIMOS pipeline (Newman et al., in preparation) esti- 
mate for a given pixel and redoing all centroid measure- 
ments. 
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Fig. 3. — Segments of the zCOSMOS (top, green), SDSS (middle, orange), and DEIMOS (bottom, black) spectra and the corresponding 
BPT diagram for COSMOS J100043. 15+020637.2. Each spectrum is shifted to the rest frame of the host galaxy, and the dashed vertical lines 
show the expected wavelengths of H/3 and [O III] A5007 (left panel) and Ha and [N II] A6584 (middle panel). For clarity, the spectra are 
offset from each other vertically and normalized to span 300 counts hr _1 . The plotted DEIMOS spectrum has been smoothed by taking the 
inverse- variance- weighted mean in a rolling 4 A wide window. The right panel shows the BPT diagram of the line flux ratios for zCOSMOS 
(green circle), SDSS (orange diamond), and DEIMOS redward (red square) and blueward (blue square) emission com ponents. The dashe d 
curve illustrates the theoretical maximum for starbursts, based on the upper limit of stellar photoionization models (Kcwlcy ct al. 2001), 
while the solid curve shows the empirical divisio n between galaxies tha t are purely star-forming and "composite" galaxies whose spectra 
are dominated by both star formation and AGN (Kauffmarm ct al. 2003). Pure star- forming galaxies lie under the dashed curve, composite 
galaxies lie between the dashed and solid curves, and pure AGN galaxies lie above the solid curve. Our line flux ratio measurements of 
COSMOS J100043.15+020637.2 show that both of its emission components are dominated exclusively by AGN activity. 



We find the two AGN emission components have pro- 
jected separations of 1.5 ± 0.5 h^ 1 kpc, 0.97 ± 0.42 h~ x 
kpc, and 1.6 ± 0.4 /i" 1 kpc in [O III] A5007, Ha, and 
[N II] A6584, respectively. The low spatial separation 
measured for the Ha emission is due to the imperfectly- 
subtracted night sky line partially obscuring the blue- 
ward portion of the Ha emission, but the [O III] A5007 
and Ha spatial offsets are roughly consistent with the 
spatial offset measured in the HST image (§ 12. ip . We 
note that the 3°. 7 difference between the position angle 
of the DEIMOS slit and the orientation of the two AGN 
(measured in § l2.1l from the HST image) produces a neg- 
ligible (0.2%) difference between the DEIMOS and HST 
measurements of projected spatial offsets. 

2.3. One- Dimensional zCOSMOS, SDSS, and 
DEIMOS Spectra 

We analyze spectra of COSMOS J100043. 15+020637.2 
from the zCOS MOS spectroscopi c redshift survey of the 
COSMOS field (jLillv et al J 120071 ) . SDSS, and DEIMOS 
to both determine the galaxy's redshift and measure 
emission line flux ratios to determine the source of its 
line emission. 

To measure the redshift of COSMOS 
J100043. 15+020637.2, we determine the value of z 
which minimizes \ 2 when comparing to a template 
spectrum. We mask out all emission lines, then fit a con- 
tinuum template spectrum based on lBruzual fc Charlotl 
(2003) stellar-population sy nthesis mo d els. The tem- 
plate spectrum, described in lYan et all (|2006T> . consists 
of a 0.3 Gyr, solar metallicity, young stellar population 
combined with a 7 Gyr, solar metallicity, old stellar 
population. From the template spectrum fits, we find 
that the zCOSMOS, SDSS, and DEIMOS spectra all 
give consistent redshifts of z = 0.36. 

The source of line emission in a galaxy is commonly 



identified using the Bal dwin-Phillips- Terlevich (BPT) di- 
agram of line ratios (|Baldwin et al.l 119811 : iKewlev et al.l 
2006). To determine the source of the line emission 
in COSMOS J100043.15+020637.2, we examine its H/3, 
[O III] A5007, Ha, and [N II] A6584 emission lines in the 
zCOSMOS, SDSS, and DEIMOS spectra (Figure^. 

The higher spectral resolution (R ~ 3000) of the 
DEIMOS spectrum enables us to discern substructure in 
the emission lines that is unresolved in the zCOSMOS 
and SDSS spectra (R ~ 600 and R ~ 1800, respec- 
tively). To determine the velocity difference between 
the two peaks we fit two Gaussians to the continuum- 
subtracted [O III] A5007 line profile, which is the emis- 
sion line with the highest signal-to-noise ratio. Based on 
the wavelengths of the peaks of the best-fit Gaussians, 
we find the line-of-sight velocity difference between the 
double peaks is 150 ± 40 km s _1 . The error in velocity 
is derived from the errors in the peak wavelengths of the 
best-fit Gaussians added in quadrature. 

To identify whether the double-peaked lines corre- 
spond to two AGN, we examine the line flux ratios 
of each emission component separately. We fit two 
Gaussians to each double-peaked line, setting the ve- 
locity separation of the peaks to the best-fit veloc- 
ity difference found for the [O III] A5007 line profile, 
but allowing the heights and widths of the Gaussians 
to vary. The areas under the best-fit Gaussians pro- 
vide estimates of the line fluxes for each emission com- 
ponent, and we find [O III] A5007/H/3=4.4 ± 1.2 and 
[N II] A6584/Ha=0.82 ± 0.17 for the blueward emis- 
sion component and [O III] A5007/H/3=3.4 ± 1.0 and 
[N II] A6584/Ha=0.85 ± 0.19 for the redward emission 
component, where the uncertainties are derived from pro- 
pogation of errors in the parameters of the best-fit Gaus- 
sians. As Figure [3] shows, the locations of these line flux 
ratios on the BPT diagram indicate that the line emis- 
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sion for both emission components is cle arly produced 
by A GN activity and not star formation (|Kewlev et al.l 
I2006H . In other words, the combined DEIMOS line 
flux ratios indicate that COSMOS J100043.15+020637.2 
hosts two distinct AGN at nearly the same redshifts. 

To measure the line flux ratios for the lower-resolution 
zCOSMOS and SDSS spectra, we subtract the contin- 
uum from each spectrum and then measure the flux of 
the H/3, [O III] A5007, Ha, and [N II] A6584 emission 
lines. For zCOSMOS we measure [O III] A5007/H/?=4.1 
and [N II] A6584/Ha=0.42 (we cannot compute errors 
on these ratios because no error array is provided 
with public zCOSMOS spectra; however, the error on 
[N II] A6584/Ha may be large because of the night 
sky line partially obscuring the Ha emission), and 
for SDSS we measure [O III] A5007/H/?=3.6 ± 0.7 and 
[Nil] A6584/Ha=0.63±0.06. The locations of these line 
flux ratios on the BPT diagram both confirm that the 
line emission in COSMOS J100043. 15+020637.2 is pro- 
duced by AGN activity (Figure [3]). 

2.4. Multiwavelength Detections 

We use the Infrared Science Archive 1 to obtain ad- 
ditional observations of COSMOS J100043.15+020637.2 
across different wavebands, which bolster the evidence 
for AGN activity found in § 12.11 and § 12.31 as described 
below. 

XMM-Newton detected fluxes of 3.82 x 10~ 14 erg cm" 2 



6.08 x 10 



-14 —2 

erg cm s 



-14 



and 2.87 x KT" erg 
cm" 2 s" 1 at energy bands 0.5 - 2.0 keV, 2.0 - 10.0 
keV, and 5.0 - 10.0 keV, respectively. These fluxes are 
above the limit for AGN detection, given by ~ 10 -15 erg 
cm" 2 s " 1 in the 0.5 - 2 keV or 2.0 - 10.0 keV energy 
bands (|Brusa et alJl2007tlCappelluti et al]l2007| ). In ad- 
dition, the Very Large Array (VLA) detected a flux of 
0.139 mJy at an observing frequency of 1.4 GHz, which is 
above the limit of 0.1 m Jy at 1.4 GHz for AGN detection 
(Schi nnerer et a l. 2007). Finally, the Spitzer Infrared Ar- 
ray Camera (IRAC) measured Vega magnitudes of 14.92, 
14.34, 13.7 9, and 12.66 at 3.6, 4.5, 5.8, and 8.0 fan, re- 
spectively ([Sanders et al.l[2007[ ). and these measurements 
fall wi thin the observed range of mid-infrared colors of 
AGN dStern et al.ll2005l ). 

These observations, combined with a dditional data 
frorn Spitzer MIPS (at 24 and 70 fim; iSanders etai] 
2007), the Canada- France-Hawaii Telescope Legacy 
Survey 2 (CFHTLS; in u*, g' , r', i', and z% SDSS (in 
u', g', r' , i' , and z'), and the Galaxy Evolution Explorer 
(GALEX; in the near-ultraviolet and far-ultraviolet) en- 
able us to compile a rest-frame SED for COSMOS 
J100043.15+020637.2 across several wavebands (Fig- 
ure [3} . This figure also shows the me dian radio-loud an d 
radio- quiet QSO SEDs compile d by lElvis et~aT1 (|1994h . 
and a lBruzual fe Chariot! (|2003l ) template SED for a 100- 
Myr-old simple stellar population with metallicity of 0.02 
times the Solar value. The QSO templates have been nor- 
malized to the 24^m detection and the stellar template 
has been normalized to the z'-band detection. The SED 
of this object is broadly consistent with a QSO that is 
obscured by dust in the optical and UV, allowing stellar 
emission from the host galaxy to dominate, with strong 

1 http://irsa.ipac.caltech.edu/ 

2 http://www.cfht.hawaii.edu/Science/CFHLS/ 
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Fig. 4.— SED of COSMOS J100043.15+020637.2, from radio 
to X-ray. Data points show fluxes (or upper limits) measured 
by various instruments, as noted in the legend. All frequencies 
hav e been shifted to the host galaxy's rest frame. Also shown are 
the Elvis et al. ( 1994) median SEDs for radio-loud and radio-quiet 
QSOs (blue and purple, respectively) nor malized to the 24-micron 
flux, as well as a lBruzual fc Charlotl ifetjul ) template SED for a 100- 
Myr-old stellar population (red) normalized to the z'-band flux. 

re-radiation by dust in the mid-infrared. Some or all of 
the infrared radiation might also arise from ongoing star 
formation, as a merger like this one would be expected to 
produce strong nuclear starbursts, but the X-ray and ra- 
dio detections confirm the presence of at least one active 
nucleus in this system. 

3. CONCLUSIONS 

The observations discussed here provide strong evi- 
dence for a dual AGN with 1.75 ± 0.03 hr 1 kpc pro- 
jected spatial separation and 150 ± 40 km s _1 line-of- 
sight velocity separation in the z = 0.36 galaxy COS- 
MOS J100043. 15+020637.2. The [O III] A5007/H/3 and 
[N II] A6584/Ha line flux ratios of the DEIMOS double- 
peaked emission lines are solid indicators that the galaxy 
hosts two AGN, a conclusion that is bolstered by the 
AGN-like luminosities of the galaxy's two central point 
sources. Line flux ratios of the unresolved, combined 
emission lines in zCOSMOS and SDSS spectra, as well 
as X-ray, radio, and infrared detections, provide further 
confirmation of AGN activity in the galaxy. Finally, the 
tidal tail visible in HST imaging is an unmistakable sig- 
nature of a galaxy merger. We conclude that COSMOS 
J100043. 15+020637. 2 is a merger-remnant galaxy with 
two inspiralling supermassive black holes, each of which 
powers an AGN. 

The discovery of this dual AGN adds significantly to 
the number of such known objects. A statistical sam- 
ple of dual AGN would provide a direct observational 
probe of both the galaxy merger rate and the kinematics 
of SMBH mergers, which are expected to produce grav- 
ity wa ves observable by n ext-generation projects such as 
LISA (|Bender et alJll998h . 

We initially identified COSMOS J100043. 15+020637.2 
as a dual AGN candidate because of its two bright cen- 
tral sources visible in HST imaging. Our finding was 
serendipitous, and there are likely more dual AGN to 
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be discovered in COSMOS. We have demonstrated that 
dual AGN candidates can be selected as bright double 
sources in EST imaging and confirmed through optical 
spectroscopy and multiwavelength observations. 
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